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Abstract 

The  results  of  a  computational  and  experimental 
study  of  feedback  control  of  Karman  vortex  shedding 
behind  a  circular  cylinder  are  presented.  The  experi¬ 
ments  were  performed  in  a  nominally  two  dimensional 
flow  around  a  12.7mm  cylinder  at  Reynolds  numbers 
up  to  10,000.  The  computations  utilized  the  two- 
dimensional  Navier-Stokes  equations  at  Reynolds  num¬ 
bers  up  to  100.  In  both  experimental  and  numerical 
studies  blowing  and  suction  slots  were  used  to  con¬ 
trol  the  vortex  shedding.  In  the  numerical  study  com¬ 
plete  suppression  of  vortex  shedding  is  achieved  up  to  a 
Reynolds  number  Re  =  80  (about  70%  above  the  onset 
Reynolds  number  for  vortex  shedding).  In  the  exper¬ 
imental  study  using  the  fluctuating  lift  force  as  feed¬ 
back,  about  50%  suppression  of  the  rms  fluctuations 
h^ls  been  achieved  at  a  Reynolds  number  of  approxi¬ 
mately  10,000. 

Nomenclature 

Cl  lift  coefficient,  L/{Q.5pdlUoo^) 

d  cylinder  diameter 

ei  unit  vector  in  the  x  direction 

ej  unit  vector  in  the  y  direction 

/  frequency  (Hz) 

ib(t)  feedback  function 

L  lift  force 

I  cylinder  span 

n  vortex  shedding  frequency  (Hz) 

Re  Reynolds  number,  Uood/v 

St  Strouhal  number,  nd/Uoo 

t  time 

u  velocity  vector  normalized  by  (/„ 

Ur  blowing/suction  normalized  velocity  (Numerical 
solution) 

Uoa  freestream  velocity 
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V 

peak  velocity  of  blowing/suction 

(experiment) 

coordinate  system,  centered  on 
aligned  with  Uoo 

cylinder  axis,  x 

Xs 

vertical  velocity  sensor  location 

a 

feedback  gain 

1/ 

kinematic  viscosity 

P 

fluid  density 

$ 

aizimuthal  angle 

T 

time  (prior  to  t) 

1  Introduction 

Many  flows  of  engineering  interest  (e.g.  separating 
flows)  produce  the  phenomena  of  vortex  shedding  and 
the  associated  structural  response.  Applications  in¬ 
clude  marine  structures,  underwater  acoustics,  civil  and 
wind  engineering  to  name  a  few.  The  ability  to  tailor 
the  wake  of  a  bluff  body  could  be  used  to  reduce  drag, 
increase  mixing  or  heat  transfer  and  enhance  combus¬ 
tion.  Vortex  shedding  from  bluff  bodies  has  been  re¬ 
viewed  most  recently  by  Sarpkaya  [1]  and  Bearman  [2]. 
Vortex  shedding  from  a  cylinder  occurs  over  wide  range 
of  Reynolds  number  from  47  <  i?e  <  300,000  and 
then  again  at  Re  >  3x10®.  From  a  mechanistic  (per¬ 
haps  simplistic)  point  of  view,  a  flow-induced  feedback 
causes  the  shedding  of  the  vortices,  where  the  chang¬ 
ing  location  of  boundary  layer  separation  on  one  side 
of  the  cylinder  causes  a  change  in  pressure  distribution 
that  forces  an  opposite  asymmetric  change  in  bound¬ 
ary  separation.  Any  small  flow  perturbation  is  enough 
to  cause  the  initial  asymmetry  that  starts  the  feedback 
mechanism.  The  observation  that  the  vortex  shedding 
is  caused  by  the  moving  location  of  boundary  layer  sep¬ 
aration  is  supported  by  the  absence  of  a  distinct  shed¬ 
ding  frequency  in  the  range  300,000  <  Re  <  3x10®.  In 
this  regime,  apparently  laminar  separation  followed  by 
turbulent  reattachment  fixes  the  turbulent  separation 
location  so  firmly  that  feedback  is  inhibited. 

The  important  aspect  of  this  phenomenon,  from  a 
control  standpoint,  is  that  the  large  pressure  forces, 
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both  lift  and  drag,  are  a  result  of  the  small  viscous 
forces  acting  in  the  boundary  layer  of  the  cylinder.  This 
implies  that  small  control  forces  applied  to  the  bound¬ 
ary  layer  can  be  used  to  control  the  larger  pressure 
forces  to  achieve  the  effect  of  modifying  the  vortex  gen¬ 
eration  and  unsteady  lift  forces. 

There  have  been  many  papers  on  the  subject  of  ac¬ 
tive.  open  loop  control  of  vortex  shedding  behind  a  bluff 
body  such  as  Koopmann  [3],  Griffin  &  Ramberg  [4], 
Armstrong  et  al.  [5]  and  Ongoren  tz.  Rockwell  [6,  7] 
to  name  a  few.  Vortex  shedding  lock-in  usually  oc¬ 
curs  when  a  cylinder  is  oscillated  in-line  to  the  flow 
direction  near  the  natural  vortex  shedding  frequency. 
However,  for  the  transverse  oscillation  of  the  cylinder, 
vortex  shedding  lock-in  occurs  near  twice  the  natural 
vortex  shedding  frequency.  Vortex  lock-in  and  con¬ 
trol  of  the  near-wake  can  also  be  observed  with  small 
rotational  oscillation  of  a  circular  cylinder  as  demon¬ 
strated  by  Tokumaru  ic  Dimotakis  [8],  Filler  et  al.  [9], 
Mo  [10]  and  Chen  et  al.  [11].  Blevins  [12],  Detemple- 
Laake  it  Eckelmann  [13],  Ffowcs  Williams  &  Zhao  [14] 
and  Roussopoulos  [15]  have  used  acoustic  forcing  as  a 
means  to  control  vortex  shedding. 

While  there  has  been  an  enormous  amount  of  effort 
put  into  the  study  of  the  open  loop  control  of  vortex 
shedding,  very  little  research  has  been  done  in  the  area 
of  the  feedback  control  of  vortex  shedding.  Monke- 
witz  [16]  investigated  the  possibility  for  the  feedback 
control  of  global  oscillations  using  a  Ginzburg-Landau 
model  equation  that  mimics  the  fluid  system.  He  finds  a 
“gain  window”  where  a  feedback  gain  with  an  optimum 
phase  is  sufficient  to  suppress  the  growth  of  the  most 
unstable  mode,  but  is  still  below  the  critical  value  at 
which  higher  modes  are  destabilized.  This  window  nar¬ 
rows  as  the  degree  of  the  supercriticality  increases  and 
control  is  typically  lost  not  far  beyond  critical  condi¬ 
tions.  Therefore,  at  higher  supercritical  Reynolds  num¬ 
bers,  the  feedback  necessary  to  stabilize  the  most  un¬ 
stable  mode  might  destabilize  the  next  unstable  mode. 
Monkewitz  et  al.  [17]  suppressed  vortex  shedding  at 
Reynolds  number  only  about  5%  above  the  onset  of  vor¬ 
tex  shedding.  Recently,  Roussopoulos  [15]  conducted  a 
similar  feedback  control  experiment  but  was  only  able 
to  suppress  the  vortex  shedding  at  Reynolds  numbers 
about  20%  above  the  onset. 

This  paper  describes  two  separate  approaches  to  in¬ 
vestigating  active  feedback  control  of  cylinder  wakes. 
The  first  approach  utilizes  the  computation  of  a  feed¬ 
back  control  of  Karman  vortex  shedding  behind  a  cir¬ 
cular  cylinder  using  a  single  downstream  sensor  and  a 
pair  of  blowing/suction  actuators  at  ±110®  from  the 
forward  stagnation  line.  The  second  effort  is  an  exper¬ 
iment  using  feedback  from  a  load  cell  which  measures 
the  oscillating  lift  forces  on  a  cylinder  to  control  a  pair 


of  blowing/suction  slots  (at  ±67  .5°)  in  flow  of  a  nomi¬ 
nal  Reynolds  number  of  10,000. 

2  Numerical  Formulation 

We  solve  the  following  Navier-Stokes  equations: 

du  2 

+  u  •  Vu  =  -Vp  + -— V-u,  (1) 

at  lie 

and  the  continuity  equation: 

Vu  =  0,  (2) 

where  u  is  the  normalized  velocity  vector,  and  p  is  the 
normalized  pressure. 

The  outer  boundary  condition  is 

u  =  e\  as  |a;|  — <■  oo,  (3) 

and  the  blowing/suction  boundary  condition  at  the 
cylinder  surface  is 

u|an  =  k(t)ur(^)(sin  fle"!  +  cos^ej),  (4) 

where  e\  and  ej  are  the  unit  vectors  in  the  x-direction 
and  y-direction,  respectively.  The  constraint  on  ur(^) 


where  0  is  the  azimuthal  angle  (i.e.  mass  is  conserved). 
We  chose  the  antisymmetrical  blowing/suction  veloc¬ 
ity  profile  with  respect  to  the  horizontal  centerline  to 
control  the  vortex  shedding  mode  which  is  antisymmet¬ 
rical. 

The  temporal  accuracy  of  the  scheme  is  second  or¬ 
der.  The  computational  grid  is  uniformly  stretched  in 
the  radial  direction  such  that  the  grid  is  denser  near 
the  cylinder  surface  to  resolve  the  high  gradient  of  the 
flow  near  the  cylinder  surface.  The  outer  computational 
boundary  is  located  at  roo  =  116.  The  number  of  grid 
points  in  circumferential  and  radial  directions  are  256 
and  150,  respectively.  The  velocities  and  pressure  are 
expanded  in  finite  Fourier  series  in  the  circumferential 
direction  and  a  second  order  accurate  finite  difference 
method  is  used  in  the  radial  direction.  Some  of  the 
solutions  in  section  4  are  presented  in  Park  et  al.  [18]. 
The  readers  are  referred  to  this  reference  and  Karni- 
adakis  et  al.  [19]  for  further  details  of  the  numerical 
method. 

3  Experimental  Procedure 

The  experiments  for  this  study  were  performed  in  the 
NCCOSC  12  inch  (305mm)  Water  Tunnel.  This  tunnel 
has  a  round  open  jet  inside  an  octagonal  test  section. 
Models  of  up  to  0.8m  length  can  be  accommodated, 
with  flow  velocities  up  to  15m/sec.  The  tunnel  also  has 
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a  24inch  (0.6m)  butterfly  valve  which  is  installed  just 
before  the  pump.  In  this  study  the  butterfly  valve  was 
used  to  throttle  the  flow  so  that  lest  section  velocities 
as  low  as  O.lm/sec  could  be  achieved. 

The  oscillatory  lift  force  on  the  test  cylinder  was  cho¬ 
sen  as  the  feedback  for  the  control  system.  A  water  tun¬ 
nel  is  the  ideal  flow  facility  for  this  study  since  the  high 
density  of  water  produces  relatively  large  lift  forces  at 
manageably  low  shedding  frequencies.  Figure  1  shows 
the  test  cylinder  apparatus  installed  at  the  exit  of  the 
open  jet  in  the  NCCOSC  water  tunnel.  The  cylinder 
was  12.7mm  in  diameter.  Lift  forces  were  measured  at 
each  end  of  the  cylinder  (top  and  bottom)  with  cus¬ 
tom  fabricated  load  beams.  Each  beam  formed  a  flex¬ 
ure  such  that  bending  the  beam  in  an  “S”  shape  would 
strain  the  opposite  legs  of  a  four  leg  Wheatstone  bridge 
strain  gage  mounted  on  the  beam.  Each  load  beam  was 
assembled  in  a  water  proof  housing,  with  the  active  end 
of  the  beam  mounted  to  the  end  of  the  cylinder  through 
a  rubber  rolling  seal.  The  rolling  seal  acted  to  keep  wa¬ 
ter  out  of  the  housing  while  permitting  a  small  amount 
of  motion  to  deflect  the  load  beam.  At  full  load  the 
deflection  of  the  cylinder  (perpendicular  to  the  flow  di¬ 
rection)  was  measured  to  be  less  than  0.04mm  or  about 
0.3%  of  the  cylinder  diameter.  The  results  presented  in 
section  5  were  at  flow  speeds  such  that  lift  loads  were 
about  0.1  of  the  maximum  load,  so  that  in  this  case, 
the  movement  of  the  cylinder  in  response  to  the  flow  is 
so  slight  as  to  be  considered  insignificant. 

Flow  visualization,  using  both  hydrogen  bubble  wire 
and  fluorescein  dye  with  Argon  laser  (514nm)  illumina¬ 
tion,  showed  that  the  flow  was  largely  coherent  across 
the  span,  however  occasionally  a  “dislocation”  in  the 
shed  vortex  could  be  seen.  Here  “dislocation”  refers 
to  a  spanwise  phase  discontinuity  of  the  shed  vortex. 
Later  tests  using  the  load  cells  showed  a  high  coherence 
between  the  lift  forces  measured  at  the  top  and  bot¬ 
tom  load  cells  (approximately  80%).  The  aspect  ratio 
(span/diameter)  of  the  test  cylinder  was  22.  With  the 
pictured  geometry  the  Strouhal  number  of  the  vortex 
shedding  from  the  test  cylinder  was  0.197  and  was  es¬ 
sentially  independent  of  Reynolds  number  in  the  range 
1,000-10,000.  This  is  close  to  but  slightly  lower  than 
the  expected  value  (0.205).  This  difference  was  prob¬ 
ably  due  to  the  blockage  effect  of  the  cylinder  in  the 
round  jet  and  the  tendency  of  the  Strouhal  number  to 
decrease  with  lower  cylinder  aspect  ratio  (Karniadakis 
&  Triantafyllou)[20]. 

Suction/blowing  slots  were  chosen  as  actuators.  Ro¬ 
tary  oscillation  may  be  more  effective  for  a  circular 
cylinder  geometry.  However,  practical  design  consid¬ 
erations  made  the  rotary  actuator  too  complex,  given 
that  lift  force  needed  to  be  measured  also.  It  was  felt 
that  the  most  effective  location  for  the  suction/blowing 


slots  was  near  to,  but  upslrearu  of  the  mean  flow  separa¬ 
tion  location  on  the  cylinder.  Numerical  icsulto  f^r  i-iie 
two-dimensional  simulation  generally  showed  this  to  be 
true.  The  mean  separation  location  is  Reynolds  num¬ 
ber  dependent,  and  for  the  low  Reynolds  number  (80) 
of  the  numerical  simulation  a  slot  location  ±110®  from 
the  stagnation  line  was  most  effective.  In  the  experi¬ 
ment,  higher  flow  velocity  corresponding  to  Reynolds 
number  in  the  9,000-10,000  range,  was  used  to  get  a 
retisonably  large  signal  on  the  load  cells.  The  low  end 
response  of  the  load  cells  was  generally  limited  by  vi¬ 
bration  from  the  water  tunnel  prime  mover  showing  up 
on  the  load  cell  output.  At  these  higher  Reynolds  num¬ 
bers  the  mean  separation  line  wanders  about  80°  from 
the  forward  stagnation  line.  A  slot  location  of  ±67.5® 
was  chosen  for  the  experimental  cylinder. 

Figure  2  is  a  cross  section  of  the  stainless  steel  cylin¬ 
der.  The  cylinder  was  constructed  in  two  pie  shaped 
sections  and  then  bolted  together.  Bolt  heads  (on  the 
aft  stagnation  line)  were  filled  and  faired  in  to  minimize 
flow  disturbance.  A  distribution  manifold  for  each  slot 
was  machined  down  the  length  of  the  larger  part,  so 
that  pressure  drop  across  the  span  of  the  slot  was  mini¬ 
mized.  Also,  flow  passages  feed  in  from  both  ends  of  the 
cylinder  minimizing  spanwise  flow  variations.  The  slots 
were  created  by  undercutting  the  smaller  part  0.076mm 
before  assembling  the  two  halves.  Imperfect  machining 
created  a  slot  width  variation  along  the  span  of  the 
cylinder  from  0.05mm  to  0.10mm,  with  the  worst  of 
the  variation  at  the  ends  of  the  cylinder. 

Suction/blowing  was  done  by  plumbing  the  slot  man¬ 
ifolds  to  a  peristaltic  pump  located  outside  the  tunnel 
test  section.  The  pump  was  connected  to  a  DC  servo¬ 
motor  such  that  a  sine  wave  input  to  the  servoamplifier 
would  produce  an  oscillatory  rotation  of  the  pump.  The 
servomotor/amplifier  acted  as  a  velocity  servo  such  that 
a  commanded  voltage  would  produce  a  corresponding 
rotational  velocity  of  the  motor.  The  pump  motor  com¬ 
bination  had  a  generally  direct  proportional  response 
up  to  30Hz.  The  plumbing  was  arranged  in  such  a  way 
that  the  amount  of  fluid  suctioned  from  one  slot  would 
be  blown  out  the  opposite  slot.  Plumbing  consisted  of 
as  much  stiff  wall  tubing  as  possible,  to  minimize  the 
capacitive  effect  of  flexible  wall  tubing.  The  peristaltic 
pump  had  three  pinch  rollers  bearing  on  the  vinyl  tub¬ 
ing  that  formed  the  flow  circuit.  Some  flow  variation 
would  occur  depending  on  the  location  of  the  pinch 
rollers  on  the  tube  relative  to  the  inlet/outiet  portion 
of  the  pump.  This  variation  was  estimated  to  be  as 
much  as  25%  for  the  open  loop  tests  presented  in  section 
5.  This  variation  was  averaged  over  20  spectrums  per 
forcing  frequency  for  the  open  loop  results  presented 
in  figures  6  through  10.  For  the  closed  loop  tests,  a 
position  transducer  was  added  to  the  pump  shaft  so 
that  the  mean  location  of  the  pinch  roller  could  be  con- 
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Figure  1:  Side  and  end  view  of  experimental  cylinder  in  NCCOSC  Water  Tbnnel.  Dimensions  in  mm. 

trolled.  Generally  the  commanded  velocity  resulted  in 
only  a  few  degrees  (approx.  10)  of  shaft  rotation  so  the 
commanded  flow  is  generally  directly  proportional  to 
shaft  rotation  in  the  closed  loop  case. 

4  Numerical  Results 

Accuracy  of  the  numerical  method  was  validated  for 
the  natural  vortex  shedding  case  by  comparing  the 
numerically  obtained  Strouhal  number  St  to  experi¬ 
mental  values  at  Reynolds  numbers  60,  100,  150  and 
200.  Agreement  between  our  results  and  the  universal 


and  v{x  =  x,,y  =  0,()  is  the  vertical  velocity  mea¬ 
sured  at  x=ij,  t/=0  at  time  t,and  a  is  the  proportional 
feedback  gain.  Changing  the  sensor  location  x,  corre¬ 
sponds  to  changing  the  phase  of  the  feedback  signal. 
The  desired  state  or  goal  is  to  drive  or  control  the  flow 
to  the  symmetric  state  with  respect  to  the  horizontal 
centerline,  which  is  the  Fornberg’s  [22]  steady  solution. 
Therefore,  we  chose  this  forcing  function  using  the  ver¬ 
tical  velocity  at  the  horizontal  centerline  since  as  the 
desired  state  (symmetric  state)  is  approached,  the  re¬ 
quired  forcing  becomes  small. 

The  location  of  the  actuators  or  the  blowing/suction 
slots  is  a  critical  factor  in  their  effectiveness  for  the 
suppression  of  the  vortex  shedding.  Blowing/suction 
slots  at  ±50",  ±65", and  ±110"  from  the  front  stagna¬ 
tion  of  the  cylinder  were  tested.  With  the  actuators 


St  —  Re  curve  of  Williamson  [21]  is  excellent  (within 

2%). 

The  feedback  function  k(t)  is  given  by  the  following 
formulation; 


k{i)  =  Q 


u(z  =  T,,y  =  0,t) 

Umor  (t) 


where 


Vmor(t)  =  maLx|v(z  =  r,,y  =  0,t)|,  (7) 

r<t 
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at  ±50°,  or  ±65”,  suppression  of  the  vortex  shedding 
rias  not  obtained.  We  postulate  the  actuators  were  in¬ 
effective  because  they  were  located  too  far  from  the 
separation  point  on  the  cylinder  surface.  The  separa¬ 
tion  points  on  the  cylinder  at  low  Reynolds  numbers 
wander  around  ±120®  as  the  vortices  are  shed.  When 
the  blowing/suction  slots  were  located  at  ±110*  from 
the  stagnation  point,  that  is  slightly  ahead  of  the  sep¬ 
aration  points  of  the  unforced  wake,  they  proved  to  be 
strong  and  effective  actuators. 

We  have  achieved  the  complete  suppression  of  vor¬ 
tex  shedding  using  feedback  up  to  Reynolds  number 
/Zc=80.  The  critical  Reynolds  number  for  the  onset  of 
the  vortex  shedding  is  47.  Reynolds  number  Re=8Q,  is 
about  70%  above  the  onset  Reynolds  number  and  cor¬ 
responds  to  the  strongly  nonlinear  state.  Monkewitz  et 
al.  [17]  and  Roussopoulos  [15]  have  achieved  complete 
suppression  of  vortex  shedding  only  up  to  Reynolds 
numbers  which  are  about  5%  and  10%  above  the  critical 
onset  Reynolds  number,  respectively. 

At  Reynolds  number  Re=80  the  suppression  of  the 
vortex  shedding  was  achieved  with  feedback  when  the 
sensor  was  located  at  x,=5.5  and  the  proportional  feed¬ 
back  gain  or=0.4.  The  optimum  proportional  feed¬ 
back  gain  a  also  depends  on  the  width  of  the  blow¬ 
ing/suction  slots.  The  vortex  shedding  suppression  was 
possible  in  the  narrow  range  5.5  —  <  <  *,  <  5.5  +  <, 
0.4  —  <  ot  <  0.4  ±  S,  where  c  and  6  arc  small.  On 

the  other  hand,  by  reversing  the  phase  of  the  feedback 
signal  (a=— 0.4),  keeping  x,  same,  we  were  able  to  en¬ 
hance  the  vortex  shedding  strength  and  increase  the  lift 
force  fluctuation  amplitude. 

Figure  3  is  the  temporal  plot  of  the  lift  coefficient  and 
the  feedback  forcing  signal.  The  time  t=0  corresponds 
to  the  fully  developed  natural  vortex  shedding  state. 
The  stabilizing  (vortex  shedding  suppression)  feedback 
forcing  (x,=5.5,  o=0.4)  is  turned  on  at  t=63.7.  As 
one  can  notice,  the  lift  force  fluctuation  amplitude  de¬ 
creases  from  0.23  to  0.01.  The  feedback  signal  starts 
out  at  an  amplitude  of  0.4  (blowing/suction  velocity 
at  40%  of  the  freestream  velocity),  but  goes  down  to 
about  0.005  (blowing/suction  velocity  at  0.5%  of  the 
freestream  velocity).  What  is  amazing  is  that  such  a 
small  amplitude  lift  force  fluctuation  level  can  be  main¬ 
tained  with  the  feedback  blowing/suction  at  mere  0.5% 
of  the  freestream  velocity.  As  soon  as  the  feedback 
forcing  is  turned  off  at  t=260,  the  lift  force  fluctuation 
begins  to  increase  exponentially  and  saturate  at  the 
natural  vortex  shedding  level. 

The  flow  structures  corresponding  to  vortex  shed¬ 
ding,  vortex  shedding  suppression,  and  vortex  shedding 
enhancement  are  presented  in  Figures  4a  through  4c. 
The  vorticity  plots  of  Figure  4a-c  are  all  at  the  same 


time 

Figure  3:  Lift  coefficient  at  Re=80.  Stabilizing  feed¬ 
back  forcing  (xj=5.5,  q:=0.4)  is  turned  on  at  t=63.7 
and  turned  off  at  t=260.  Destabilizing  feedback  forcing 
(x,=5.5,  a=-0.4)  is  turned  on  at  t=475.7.  a)  t=43.04, 
b)  t=170.9,  c)  t=243.58,  d)  t=452.80,  e)  t=640.3. 

phase  of  the  vortex  shedding  or  at  the  peak  of  the 
lift  force  fluctuation.  Figure  4a  is  the  natural  vortex 
shedding.  Figure  4b  is  the  reduced  voriex  shedding 
state  some  time  after  the  suppressing  feedback  forcing 
is  turned  on.  Note  that  the  vortices  being  shed  are  elon¬ 
gated  in  the  streamwise  direction.  Finally,  figure  4c  is 
the  fully  suppressed  state,  where  the  vortex  structure 
looks  almost  symmetric  top  and  bottom. 

The  results  at  slightly  higher  Reynolds  number 
Re=l00  are  significantly  different.  Figure  5  is  the  lift 
coefficient  plot  for  the  feedback  sensor  location  Xj=5.05 
and  0=0.3  with  slot  width  slightly  larger  than  that  of 
the  Re=80  case.  The  feedback  control  is  turned  on  at 
t=422.  Lift  fluctuation  decreEises  until  t=560,  but  after 
that  it  begins  to  increase  again  and  maintains  itself  at  a 
near-equilibrium  state.  The  Strouhal  frequency  of  the 
vortex  shedding  at  this  near-equilibrium  state  is  0.1171 
as  compared  to  the  natural  vortex  shedding  Strouhal 
frequency  0.1654.  This  behavior  has  been  suggested 
by  Monkewitz  in  his  theoretical  work.  He  suggested 
that  the  wake  stabilization  might  only  be  possible  near 
the  onset  of  shedding  in  a  gain  window  between  the 
suppression  of  the  first  mode  and  the  destabilization  of 
the  higher  mode.  Recently,  Roussopoulos  [15]  has  ex¬ 
perimentally  observed  a  similar  behavior  at  Reynolds 
number  Re=52.8. 

5  Experimental  Results 

The  open  loop  response  of  the  lift  to  the  forcing  of 
the  suction/blowing  actuators  was  studied  first.  The 
open  loop  forcing  response  was  principally  studied  to 
provide  some  insight  in  the  design  of  a  closed  loop  con¬ 
trol  system.  Also  of  interest  was  an  examination  of  the 
“receptivity"  and  “lock-in”  regions  for  open  loop  wake 
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Figure  4:  Vorticity  contour  plot  at  /Zc=80  correspond¬ 
ing  to  figure  3.  a)  t=43.04,  b)  t=170.9,  c)  t=243.58,  d) 
t=452.89,  c)  t=640.3. 

forcing  discussed  by  Karniadakis  &  TViantafyllou[20]. 
In  this  context  lock-in  implies  that  the  shedding  fre¬ 
quency  adjusts  itself  so  that  it  exactly  matches  the  forc¬ 
ing  frequency.  For  large  forcing  amplitudes,  the  range 
of  frequencies  over  which  the  shedding  frequency  fol¬ 
lows  the  forcing  frequency  becomes  larger.  Receptivity 
on  the  other  hand  is  defined  as  that  region  in  the  two 
dimensional  forcing  frequency-amplitude  space,  where 
the  wake  responds  to  the  forcing  in  some  manner,  but 
does  not  actually  adjust  its  frequency  to  match  the  forc¬ 
ing  frequency.  This  region  is  larger  and  surrounds  the 
lock-in  region. 

For  both  the  open  loop  and  closed  loop  tests  a 
freestream  velocity  of  0.745m/sec  was  chosen.  At  this 
velocity  the  load  cells  produced  a  sufSciently  large 
signal  for  use  as  feedback,  and  the  natural  shed¬ 
ding  frequency  was  fairly  low  at  about  11.3Hz,  which 
was  well  within  the  proportional  response  range  of 
the  pump/motor  combination.  All  results  presented 
are  nondimensionalized  using  the  quantities  nondimen- 
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Figure  5:  Lift  coefhcient  at  /2e=100.  Feedback  forcing 
(i,=5.05,  0=0.3)  is  turned  on  at  t=422. 

sional  frequency  =  fd/Uoo,  forcing  coefficient  v/Uoo, 
lift  coefficient  Cl  and  Reynolds  number  Re. 

Figure  6  is  a  waterfall  display  of  the  nondimen- 
sional  linear  spectrum  of  the  lift  coefficient,  as  the  suc¬ 
tion/blowing  forcing  frequency  is  stepwise  (81  steps) 
increased  from  fdlUoo  =  0.135  to  0.22.  Here  the  peak 
forcing  amplitude  is  v/Ueo  =  0.04,  a  relatively  small 
amount  of  forcing.  The  forcing  frequency  shows  as 
a  diagonal  peak  in  the  waterfall  display,  starting  at 
the  lower  left  corner  and  advancing  toward  the  upper 
right.  The  natural  shedding  frequency  shows  as  a  con¬ 
tinuous  ridge  at  /d/t/oo  0.197.  At  this  low  forc¬ 
ing  level(v/l/oo  =  0.04)  the  natural  shedding  is  only 
slightly  affected.  The  only  response  appears  at  the  odd 
frequency  of  fd/Uoo  ~  0.17,  where  careful  examination 
shows  a  very  small  response  in  all  the  natural  spec- 
trums.  This  data  indicates  that  this  level  of  forcing 
is  below  the  “threshold”  postulated  by  Karniadakis  ic 
'IViantafyllou[20].  Note  that  the  momentum  flux  of  the 
flow  through  the  slots  itself  will  produce  a  force  read¬ 
ing  on  the  load  ceils.  The  relative  magnitude  of  this 
reaction  force  was  tested  by  stepping  through  the  test 
frequencies  with  no  flow  in  the  test  section.  In  general 
this  reaction  force  was  at  least  an  order  of  magnitude 
less  than  the  forces  induced  with  flow  over  the  cylinder 
and  is  essentially  independent  of  frequency.  The  sec¬ 
ond  harmonic  of  the  forcing  frequency,  which  shows  in 
the  lower  right  hand  quarter  of  the  waterfall  graphs  is 
apparently  an  artifact  of  the  reaction  forces. 

At  a  higher  forcing  amplitude  in  figure  7,  as  the  forc¬ 
ing  frequency  approaches  the  natural  frequency  it  can 
be  observed  that  at  fd/Uoo  ~  0.194,  the  shedding  fre¬ 
quency  jumps  to  a  lower  frequency  to  coincide  with  the 
forcing  frequency.  This  would  be  the  left  side  of  the 
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Figure  6:  Waterfall  plot  of  lift  coefficient  spectrum  at  a 
peak  suction/blowing  forcing  of  v/Uoa  =  0.04.  Forcing 
frequency  is  discretely  swept  from  fd/Uoo  =  0.136  to 
0.22.  Each  spectrum  is  averaged  20  times. 


lock-in  region.  Not  readily  apparent  from  figure  7  (but 
nontheless  apparent  on  individual  examination  of  the 
spectrums)  is  an  area  of  sympathetic  response  at  both 
the  forcing  frequency  and  the  natural  frequency,  where 
the  response  amplitude  for  both  peaks  is  increased. 
This  may  correspond  to  the  receptive  region  of  Kar- 
niadakis  k  Triantafyllou[20].  Their  numerical  results 
indicate  that  the  induced  response  in  this  receptive 
region  was  the  shedding  frequency  jumping  back  and 
forth  between  the  forcing  and  the  natural  frequency, 
rather  than  the  coexistence  of  two  distinct  frequencies. 
Our  observations  did  not  suggest  this,  as  the  time  trace 
of  the  lift  (in  this  response  region),  looked  very  much 
like  the  typical  superposition  of  two  frequencies.  This 
would  generally  be  a  modulated  sine  wave  with  the 
modulation  (“beat”)  frequency  equal  to  the  difference 
of  the  two  basic  frequencies.  One  further  observation 
is  that  at  forcing  frequencies  higher  than  natural,  there 
is  almost  no  response  at  all. 

Figures  8  and  9  are  similar  to  figure  7  with  yet  higher 
amplitude  forcing.  The  trends  at  the  higher  forcing 
level  are  consistent  with  previous  observations.  At 
these  higher  forcing  levels  the  lock-in  and  receptive  re¬ 
gions  are  larger.  However,  there  is  still  little  response  to 
forcing  at  higher  than  the  natural  frequency.  In  figure  9 
close  observation  of  the  natural  frequency  ridge  shows 
a  slight  bending  of  the  natural  frequency,  towards  the 
forcing  frequency,  as  the  lock-in  region  is  approached 
from  below.  Figure  10  shows  the  “frequency  response” 
of  this  system  at  two  forcing  levels  {v/Uoo  =  0.06  and 
vfUoo  =  0.14).  Here  the  term  frequency  response  is 
used  in  the  linear  system  sense.  The  curves  in  fig¬ 
ure  10  represent  the  output  of  the  system  (Cl  at  the 
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Figure  7:  Same  as  figure  6,  with  suction/blowing  forc¬ 
ing  of  v/Uoo  =  0.079. 
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Figure  8:  Same  as  figure  6,  with  suction/blowing  forc¬ 
ing  of  v/Uoo  =  0.119. 

forcing  frequency)  divided  (as  complex  numbers)  by 
the  input  (v/Uoo)  as  a  function  of  forcing  frequency. 
The  response  is  represented  as  a  complex  number  with 
both  magnitude  and  phase.  Linear  response  systems 
are  characterized  by  magnitude  and  phase  responses 
independent  of  input  level.  The  curves  in  figure  10  rep¬ 
resent  the  “vector”  average  of  20  spectrums  for  each 
of  81  response  values.  Vector  averaging  is  analogous 
to  ensemble  averaging  of  time  series,  except  performed 
in  the  frequency  domain.  Interesting  aspects  of  fig¬ 
ure  10  are:  1)  The  change  of  phase  of  the  response  is 
independent  of  forcing  level.  2)  Phase  change  is  large 
in  the  region  near  the  natural  frequency.  3)  Ampli¬ 
tude  response  is  not  directly  proportional  to  forcing 
(not  surprising  for  a  nonlinear  system).  4)  The  peak 
response  for  the  lower  level  forcing  is  close  to  the  natu- 


7 


C. 


0.10  0.20  0.30 

fd/U. 

Figure  9:  Setme  as  figure  6,  with  suction/blowing  forc¬ 
ing  of  v/Uoo  =  0.158. 


Figure  10:  Lift  coefficient  response  CiUaah  (ampli¬ 
tude  and  phase)  at  the  forcing  frequency,  with  suc¬ 
tion/blowing  peak  forcing  velocity  of  v/(7oo  =  0.14  and 
0.06. 

ral  frequency.  The  higher  level  forcing  moves  the  peak 
response  to  a  slightly  lower  frequency.  This  is  the  bend¬ 
ing  of  the  natural  frequency  referred  to  previously  in  the 
discussion  of  figure  9.  5)  A  receptive  response  region 
occurs  at  fd/Uoo  0.17. 

The  independence  of  the  phase  response  to  forc¬ 
ing  level  (noted  in  figure  10)  suggests  that  perhaps  a 
strictly  linear  feedback  law  could  be  used  to  control 
the  level  of  unsteady  lift  forces.  As  a  first  cut,  an  IBM 
clone  microcomputer  was  programmed  to  serve  as  a 
PID  (proportional,  integral  and  derivative)  controller. 
Input  to  the  control  loop  was  the  (amplified  and  an¬ 
tialiased)  average  of  the  top  and  bottom  load  cell  lift 
forces.  C  .put  was  an  analog  voltage  to  drive  the  pump 
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Figure  11:  Closed  loop  forcing  using  PD  feedback,  with 
control  turned  on  at  15  seconds  and  turned  off  at  55 
seconds. 
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Figure  12:  Response  of  lift  coefficient,  same  conditions 
as  figure  1 1 . 

servomotor.  Parameters  of  the  PID  loop  were  manually 
adjusted  to  effect  the  desired  results.  Here  the  loop  er¬ 
ror  was  taken  to  be  any  lift  force  different  from  zero. 
The  integral  term  of  the  PID  loop  was  set  to  zero,  as 
this  term  controls  the  steady  state  error  of  the  lift,  gen¬ 
erally  of  no  concern  for  this  study.  Figure  11  shows  the 
nondimensionalized  control  signal  to  the  pump,  as  the 
“best”  manually  adjusted  control  loop  is  first  turned 
on  and  then  off.  Figure  12  is  the  corresponding  long 
term  response  of  the  (top)  lift  force  as  the  the  control  is 
turned  on  and  off.  The  results  are  not  spectacular,  but 
the  reduced  level  of  unsteady  lift  in  the  middle  section 
of  the  trace  of  figure  12  represents  a  50%  reduction  of 
the  rms  level  of  lift  force.  The  parameters  of  the  con¬ 
trol  loop  also  could  be  set  to  enhance  the  unsteady  lift 
forces,  however  the  response  level  of  the  system  was  so 
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Figure  13:  Averaged  spectrum  of  lift  coefficient,  with¬ 
out  any  control.  Same  conditions  zis  the  no  control  case 
of  figure  12. 


Figure  14:  Averaged  spectrum  of  lift  coefficient,  with 
control.  Same  conditions  as  the  control  on  case  of  fig¬ 
ure  12. 


high  as  to  cause  concern  for  the  mechanical  integrity  of 
the  pump  system  and  these  results  were  not  pursued. 

Figure  13  shows  the  spectrum  of  lift  response  before 
control  is  applied  and  figure  14  is  taken  during  the  pe¬ 
riod  when  the  control  loop  is  operative.  In  figure  14  it  is 
observed  that  in  addition  to  the  reduction  of  the  level  of 
the  natural  shedding  frequency,  an  enhanced  response 
is  appearing  at  fd/Uoo  =  0.17.  In  the  process  of  the 
averaging  of  the  spectrum  of  figure  14,  it  was  observed 
that  the  system  was  alternately  jumping  between  the 
natural  frequency  and  the  new  preferred  frequency  of 
0.17.  This  is  in  contrast  to  the  open  loop  response  to 
forcing,  which  can  be  characterized  as  the  coexistence 
of  two  distinct  frequencies.  This  behavior  can  be  com¬ 
pared  with  the  numerical  results  of  figure  5.  Here  the 
flow  with  feedback  has  chosen  a  new  (lower)  shedding 
frequency.  Again,  this  behavior  is  similar  to  observed 
behavior  of  the  systems  in  Monkewitz  [16],  Monkewitz 
et  al.  [17]  and  Roussopoulos  [15].  We  believe  this  be¬ 
havior  to  be  endemic  of  linear  control  systems  applied 
to  nonlinear  systems  of  this  type,  at  high  supcrcritical- 
ity  (e.g.  high  Reynolds  number). 

6  Conclusions 

The  computational  study  of  the  feedback  control  of 
the  Karman  vortex  shedding  behind  a  circular  cylinder 
has  been  undertaken  at  low  Reynolds  numbers.  Two- 
dimensional  Navier-Stokes  equations  are  solved  numeri¬ 
cally.  The  actuators  are  a  pair  of  blowing/suction  slots 
located  at  ±110’  from  the  leading  stagnation  point. 
A  single  sensor  is  used  and  the  actuators  are  out  of 
phase  with  each  other.  Complete  suppression  of  vortex 
shedding  is  achieved  up  to  Reynolds  number  Re=80, 
which  is  about  70%  above  the  onset  Reynolds  number 


of  vortex  shedding.  At  Re=100  the  feedback  suppres¬ 
sion  destabilizes  a  second  mode. 

Closed  loop  feed  back  stabilization  has  been  used 
to  suppress  lift  fluctuation  at  Reynolds  number 
Re=  10,000  up  to  50%  over  the  natural  conditions.  Al¬ 
though  the  experiments  that  were  performed  were  done 
at  a  much  higher  Reynolds  number  than  the  simula¬ 
tions,  many  of  the  salient  features  of  the  two  studies 
were  similar.  The  open  loop  response  of  the  cylinder 
to  high  amplitude  forcing  (Figure  9)  exhibited  a  slight 
bending  of  the  natural  frequency,  towards  the  forcing 
frequency,  as  the  lock-in  region  is  approached  from  be¬ 
low.  This  observation  is  similar  to  some  numerical  re¬ 
sults  of  one  of  us  (Park)  with  open  loop  suction/blowing 
forcing  at  Re  =  120.  We  believe  this  to  be  related  to  a 
change  of  the  wake  structure  (wider),  such  that  the  pre¬ 
ferred  shedding  frequency  is  lowered.  In  addition,  the 
destabilization  of  a  second  instability  by  the  feedback 
control  was  noted  in  the  experiments  and  mirrored  in 
the  simulations. 

We  would  like  to  acknowledge  the  support  of  the 
ONT  (Office  of  Naval  Technology)  Postdoctoral  Fel¬ 
lowship  to  one  of  us  (Daniel  Park). 
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